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Abstract 

Debris dust in the habitable zones of stars - otherwise known as exozodiacal dust - 
comes from extrasolar asteroids and comets and is thus an expected part of a planetary 
system. Background flux from the Solar System’s zodiacal dust and the exozodiacal 
dust in the target system is likely to be the largest source of astrophysical noise in 
direct observations of terrestrial planets in the habitable zones of nearby stars. Fur- 
thermore, dust structures like clumps, thought to be produced by dynamical interactions 
with exoplanets, are a possible source of confusion. In this paper, we qualitatively as- 
sess the primary impact of exozodical dust on high-contrast direct imaging at optical 
wavelengths, such as would be performed with a coronagraph. Then we present the 
sensitivity of previous, current, and near-term facilities to thermal emission from debris 
dust at all distances from nearby solar-type stars, as well as our current knowledge of 
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dust levels from recent surveys. Finally, we address the other method of detecting debris 
dust, through high-contrast imaging in scattered light. This method is currently far less 
sensitive than thermal emission observations, but provides high spatial resolution for 
studying dust structures. This paper represents the first report of NASA’s Exoplanet 
Exploration Program Analysis Group (ExoPAG). 

Subject headings: Extrasolar Planets, Astronomical Techniques 


1. Introduction 

Interplanetary dust interior to the Solar System’s asteroid belt is called the zodiacal dust, 
which comes from comet comae and asteroid collisions, just like the dust in any debris disk. At 
least 16% of nearby solar-type stars harbor cold, outer debris dust much denser than the dust in 
the Solar System’s Kuiper belt (Trilling et al. 2008). Sensitive new far-IR surveys probing similar 
stars for cold dust down to about 10 times the Solar System Kuiper belt level find a detection rate 
of ~ 25% (Eiroa et al., in preparation). Unfortunately, we know little about warmer exozodiacal 
dust (exozodi) in the inner reaches of nearby systems, i.e. within the stars’ habitable zones. About 
1% of nearby solar-type stars show a large amount of emission from warm dust (> 1000 times 
the expected zodiacal dust emission in the wavelength range from 8.5 //m to 12 fim; Lawler et al. 
2009). However, a more sensitive survey for exozodiacal dust around a smaller set of nearby stars 
with the Keck Nulling Interferometer (KIN) found mostly non-detections (discussed further below; 
Millan-Gabet et al. 2011). As will be shown, background flux from our local zodiacal dust and the 
exozodi will likely dominate the signal of an Earth-analog exoplanet (an exoEarth, for short) in 
direct images and spectra, even if exozodi levels are no greater than the Solar System level. 

Furthermore, debris disks often display dust structures, such as narrow rings (e.g. HR 4796A; 
Schneider et al. 2009), clumps (e.g. e Eri; Greaves et al. 2005), and warps/inclined sub-disks (e.g. 
AU Mic; Krist et al. 2005). A notable example of a warped disk - the famous /3 Pictoris system - 
is shown in Figure 1. Such structures are thought to be produced by the dynamical influence of a 
planet, as are the clumps of zodiacal dust leading and trailing the Earth in its orbit (Dermott et al. 
1994). Dust clumps are likely to be the most troublesome source of confusion in direct imaging 
of Earth-size planets, since dynamical clumps orbit the star, though not necessarily with the same 
period as the perturbing planet (e.g. Kuchner &: Holman 2003). 

Therefore, exozodiacal dust complicates direct imaging of exoplanets in two ways: 1) as a 
source of noise and 2) as a source of confusion. The exozodical dust levels around nearby stars 
will be as important to the success of efforts to find and characterize Earth-like exoplanets as the 
fraction of stars with potentially habitable planets (r?Earth)- There is a strong need to sensitively 
probe nearby stars for small amounts of dust in their habitable zones. In this paper, we will 
first demonstrate the effect of exozodi emission as a source of increased noise in direct imaging of 
exoplanets, then assess the current knowledge of debris dust abundances and expected progress 


toward increasing our knowledge. 


1.1. Quantifying Debris Dust: What is a “Zodi?” 

Debris disk observers typically express the amount of dust in a debris disk using the system’s 
fractional dust luminosity, Tdust/D* 5 which is the light absorbed by the dust and reemitted at 
thermal wavelengths (infrared to millimeter) relative to the stellar luminosity. This parameter, 
often also called the fractional infrared luminosity, is determined by integrating the total long- 
wavelength excess flux seen in a disk spectral energy distribution (SED) over frequency. Therefore, 
it can be measured from unresolved photometry, more easily at long wavelengths where the stellar 
emission is relatively faint. L& nst /L± is not a unit of optical depth, dust mass, or surface brightness. 
In the optically thin case, it is proportional to the dust mass but is affected by grain properties like 
size and composition. 

In the context of direct exoEarth observations, the abundance of debris dust in habitable zones 
is typically given in units of “zodis,” which arose as a convenient way of quickly expressing some 
sort of ratio to the Solar System zodiacal dust. However, extreme caution must be taken when 
interpreting the unit, since different workers have used different dust properties when constructing 
their one zodi. The definition adopted largely depended on what was useful or possible for the 
analysis at hand. For example, Gaidos (1999) defined one zodi (spelled “zody”) as the effective 
emitting area of the Solar System’s zodiacal dust. Roberge & Kamp (2010) took one zodi to be the 
fractional dust luminosity of the zodiacal dust (Ld ust /L* ~ 10” 7 ; Nesvorny et al. 2010), in order to 
relate the dust abundances of known debris disks to the Solar System (by this standard, (3 Pictoris 
has about 10,000 zodis of cold dust). 

Others assume one zodi is a debris disk that is identical to the Solar System’s zodiacal dust in 
every respect, including total mass, spatial distribution, grain size distribution, albedo, scattering 
phase function, etc. (e.g. Millan-Gabet et al. 2011). Here, we will refer to such a disk as a zodiacal- 
twin disk; around a Sun-twin star, we call it a Solar System-twin disk and it has Td ust /L* ~ 10“ ' . 
However, care must be taken when applying a zodiacal-twin disk to a star different than the Sun, 
since such a disk around another type of star may not be both self-consistent and truly identical to 
the zodiacal dust in every way. ZODIPIC 1 is a publicly available code frequently used to calculate 
zodiacal-twin disks for any type of star, using the zodiacal dust model in Kelsall et al. (1998), 
which was itself based on fits to observations of the zodiacal dust emission with the COBE DIRBE 
instrument. The code does not change the dust density distribution when varying the central star, 
as expected for a zodiacal-twin disk. However, the inner disk radius set by the dust sublimation 
temperature does move outward with increasing stellar luminosity, to prevent an unphysical sit- 
uation with dust at impossibly high temperatures. The unfortunate consequence of this is that 


1 Available for download at http://asd.gsfc.nasa.gov/Marc.Kuchner/home.html. 
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the total dust mass and Ld ust /L* value of a ZODIPIC zodiacal- twin disk decreases as the stellar 
luminosity increases. 

Nonetheless, for direct observations of exoplanets at any wavelength, it is the surface bright- 
ness of the dust thermal or scattered emission near the location of the exoplanet that matters. 
While L dust /T* can be measured for debris disks, it is not straightforward to convert T du st/T* into 
a surface brightness; one must assume the spatial distribution of the dust, as well as many other 
properties that are also typically unknown (e.g. albedo). But declaring one zodi to be a zodiacal- 
twin disk is also perilous. For most debris disks, we do not know if the inner dust matches the 
zodiacal dust and in some cases we know it cannot (e.g. HD 69830; Beichman et al. 2005). Fur- 
thermore, this definition is typically applied to stars that are not identical to the Sun without any 
adjustment to the radial dust distribution, which is not likely to provide an accurate comparison to 
the Solar System (discussed further below). Because of these problems, we eschew use of the unit 
“zodi” in this paper, although we acknowledge its convenience and anticipate its continued usage. 


1.2. Exozodi Surface Brightness 

In the next section of this paper, we show the impact of exozodi emission on direct observations 
of exoEarths at optical wavelengths. To do this, we need a reference value for the scattered light 
surface brightness. For a Solar System- twin disk viewed at 60° inclination (the most likely value 
in a random distribution of inclinations), the modeled surface brightness at a projected separation 
of 1 AU in the plane of the disk (the quadrature position for a planet) and A = 0.55 gin is about 
22 mag arcsec -2 . Other brightness values calculated with ZODIPIC appear in Table l.B-1 of the 
TPF-C STDT Final Report 2 . The inclination of the system strongly affects the observed surface 
brightness; at 1 AU, the surface brightness of a ZODIPIC zodiacal-twin disk viewed edge-on is 
nearly 3 times greater than when viewed face-on. 

To apply the above surface brightness value to another star, we take advantage of the fact that 
by definition the habitable zone is the region around a star where an Earth-like planet receives the 
right amount of energy for it to have liquid water on its surface. Therefore, the habitable zone 
moves with changing stellar luminosity to keep the incident flux constant, such that 

f EEID = 1 AU * yj L*/ Lq , (1) 

where EEID stands for “Earth equivalent insolation distance.” No matter the star, bodies in the 
habitable zone always receive about the same total amount of flux; if their physical properties are 
the same, they also reflect or emit about the same total flux. 

One consequence of this is that the absolute bolometric magnitude of an Earth-twin in the 
habitable zone does not change from star to star, assuming the same viewing geometry. Another is 


http : / / exep . jpl . nasa . gov/f iles/exep/STDT_Report_Final_Ex2FF86A . pdf 
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that, assuming identical grain properties and inclination, the same amount of dust in the habitable 
zones of different stars has roughly the same surface brightness. Therefore, we can reasonably adopt 
the surface brightness at 1 AU of a Solar System-twin disk viewed at 60° inclination as a reference 
value for the exozodi brightness at the EEID for all stars. 

The movement of the habitable zone with stellar luminosity explains the aforementioned prob- 
lem with applying an unmodified zodiacal-twin disk to stars that are not Sun-like. For an earlier- 
type/later-type star, the EEID will be at a larger /smaller radius than 1 AU, probing a nrore/less 
distant region of the zodiacal-twin disk where the surface density may not be the same; in fact, the 
Kelsall et al. (1998) model says it should be lower/higher. If one was expecting the same amount 
of dust in the star’s habitable zone as in the Solar System’s, the zodiacal-twin disk would need to 
be stretched/shrunk in radius by the square-root of the stellar luminosity. 


2. Impact of Exozodical Emission on Exoplanet Direct Imaging 


The importance of zodiacal and exozodiacal background emission is shown by comparing the 
counts from both sources to the counts from an exoEarth. Note that the following discussion applies 
to ordinary high-contrast imaging, such as would be performed with a coronagraph. The impact 
of exozodi on mid-IR nulling interferometers intended for direct exoEarth observations (e.g. ESA’s 
DARWIN mission concept) is extensively discussed in Defrere et al. (2010). Following the approach 
in Brown (2005) and assuming a uniform exozodi distribution near the planet location, the counts 
ratio is 
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where n p ; x is the number of pixels in a critically sampled, diffraction-limited point-spread func- 
tion (PSF) from a partially obscured circular aperture (n p i x = 1/sharpness « 1/0.07 ~ 14.3); 
II = (jWd ^TradfaET 2 ) 2 ^ ie angular area of a pixel at wavelength A for a telescope of diameter 
D] z is the surface brightness of the local zodiacal dust, which depends on the direction to the tar- 
get star (a generic value is ~ 23 mag arcsec -2 at A = 0.55 /im); e is the exozodi surface brightness 
at the EEID relative to the brightness of a Solar System-twin disk; x is the surface brightness at 
the EEID (1 AU) of a Solar System-twin disk viewed at 60° inclination (~ 22 mag arcsec -2 for 
A = 0.55 /im); M p is the absolute magnitude of the planet ( My = 29.7 mag for an Earth-twin 
orbiting at the EEID, viewed at quadrature); and d is the distance to the system in pc. 


The left-hand panel of Figure 2 shows the counts ratio versus exozodi brightness (e), assuming 
A = 0.55 n m, d = 10 pc, the parameter values given above, and three different telescope diameters. 
For a 4-meter telescope, the counts from the local zodiacal dust plus a Solar System-twin disk of 
exozodi viewed at 60° inclination are about 5 times greater than the counts from the Earth at 10 pc. 
While this may sound disastrous, astronomers commonly observe targets that are fainter than the 
surrounding background, at the price of increased noise and longer exposure times. Assuming 
background-limited imaging and including unsuppressed light from the central star, the time to 
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detect an exoplanet at some S/N in an image is 
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where Fq is the specific flux for zero magnitude in the image bandpass (~ 9500 counts s” 1 cm” 2 nm” 1 
in V band); AA is the image bandpass width; T is the total facility throughput; A = ( 20 l^Llia iF 5 ) > 
( is the contrast level at the position of the exoplanet with respect to the theoretical Airy peak of 
the stellar image (< 10” 10 for instruments designed for direct imaging of exoEarths); and m* is the 
stellar apparent magnitude (for more details on the derivation of this equation, see Brown 2005). 


The terms within the square brackets of Equation 3 account for the background sources: the 
first for local zodiacal dust, the second for exozodi, and the third for unsuppressed starlight. The 
third term is much smaller than the sum of the other two for all known mission concepts designed 
for direct exoEarth observations at optical/near-IR wavelengths. This equation does not include 
any dark current or read noise, which are negligible compared to the other noise terms. The 
important points to take away are that in this idealized situation 1) the exposure time decreases 
as D 4 - one factor of D 2 for the collecting area and another for the smaller PSF leading to less 
background blended with the planet - and 2) the time increases linearly with increasing exozodi 
surface brightness. This is shown in the left-hand panel of Figure 2, where the imaging time is also 
plotted for the three telescope diameters. 


The goal of a mission to obtain direct exoEarths observations is to maximize the number 
of planets imaged and then characterized with spectra, while maintaining a reasonable mission 
architecture and lifetime. The right-hand panel of Figure 2 shows how ^Earth and the exozodi level 
combine to impact the performance of such a mission (ignoring confusion). The plot was created 
using a simple mission planning code that chooses real stars within 30 pc of the Sun for observation 
until the mission lifetime is reached (the mission parameters and planning code are detailed in 
Turnbull et al. 2012). It shows that the smaller ^Earth is, the lower the exozodiacal dust level that 
can be tolerated while still obtaining direct observations of a sufficient number of exoplanets over 
the mission lifetime. 


The analysis in this section demonstrates the primary effect of exozodi on direct exoEarth 
observations, i.e. increased background noise and longer exposure times. However, Equation 3 
applies exactly only in the ideal case of uniform background emission and no systematic errors (i.e. 
no speckles). Furthermore, the methods used to extract the exoplanet signal from the background 
will affect the S/N achievable. If the background contains no unresolved structures, then matched 
filtering should be an efficient way of identifying point-sources like exoplanets (see Kasdin & Braems 
2006). However, since unresolved dust clumps and other background sources (e.g. stars and galaxies) 
should often be present, this technique will not resolve all confusion issues. Another option is to 
produce a model for the exozodi emission, fitting to portions of the exoplanet image, then subtract 
the model from the image. This should remove smoothly varying exozodi emission fairly easily. To 
account for unresolved dust clumps, more complex dynamical modeling of the whole system (dust 
and exoplanets together) may be required. There are also observing strategies to mitigate sources 


of confusion (e.g. multi-epoch imaging, multi-color imaging, or even direct spectroscopy). But as 
of this writing, the problem of confusion for direct exoplanet imaging has not been thoroughly 
studied. 


3. Thermal Emission from Dust 

Now that we have established the importance of exozodi, we assess our knowledge of the 
debris dust levels around nearby solar-type stars. Assuming the stellar spectrum is well-described 
by a Rayleigh- Jeans law (as it typically is at mid-IR and longer wavelengths) and the dust thermal 
emission is single-temperature blackbody radiation, the fractional dust luminosity may be expressed 
as 

£d„„ fF dM \ 

L, ! K ) h V T ; 3 ’ ' ' 

where Edust and E* are the dust and stellar fluxes at some frequency u, T ( j is the dust temperature, 
and E* is the stellar effective temperature (e.g. Bryden et al. 2006). For blackbody grains, assuming 
all the dust is at the same temperature is equivalent to assuming that all the grains are at the same 
distance from the star (i.e. a ring- like disk). In reality, grains of different sizes will have different 
temperatures at the same distance. That being said, most debris disk SEDs are fairly well-described 
by either a cold (~ tens of K) single-temperature blackbody - an outer dust ring - or the sum of a 
cold blackbody and a warmer one - outer and inner dust rings, an architecture reminiscent of the 
Solar System’s Kuiper and asteroid belts (e.g. Chen et al. 2006, 2009). 

Using Equation 4, we may estimate the sensitivity of recent, current, and near-term telescope 
facilities to exozodiacal dust at different temperatures. Figure 3 shows the estimated 3cr sensi- 
tivity curves for the Spitzer Space Telescope (Spitzer), the Wide- field Infrared Survey Explorer 
(WISE), the Keck Interferometer Nuller (KIN), the Herschel Space Observatory (Herschel), the 
Large Binocular Telescope Interferometer (LBTI), the James Webb Space Telescope (JWST), and 
the Atacama Large Millimeter Array (ALMA). The relevant facility performance parameters - the 
instruments, years of operation (actual or expected), observation wavelengths, and uncertainties - 
appear in Table 1. 

To calculate the sensitivity curves, we assumed a stellar effective temperature equal to the 
Sun’s. For each instrument, the lcr photometric uncertainty at the observation wavelength (Col- 
umn 5 of Table 1) was converted into an uncertainty on the flux relative to the stellar flux at that 
wavelength (o"F dust /F* i n Column 6). We then replaced Ed us t/E* in Equation 4 with three times the 
relative flux uncertainties and plotted Edust/E* as a function of Td us t- Each curve has a minimum at 
the temperature for which the dust blackbody emission curve peaks at the frequency of observation 
(Tpe a k in Column 4). Emission from dust warmer/cooler than T pea k can be detected with reduced 
sensitivity. However, in most systems, emission from dust at temperatures near T pea k will swamp 
any emission from much warmer/cooler dust, which is why a survey for habitable zone dust is best 
done with an instrument that operates near 10 /im. Note that although JWST/MIRI will have 


only moderate absolute photometric accuracy, the telescope’s large collecting area will allow MIRI 
to achieve this accuracy for fainter stars than other facilities covering its bandpass. 

The temperature ranges for dust in the habitable zone and in the Kuiper belt of a Sun-like star 
are overlaid in Figure 3. The temperature range for the habitable zone is 208 - 312 K, calculated 
assuming the zone extends from 0.8 AU to 1.8 AU (Kasting et al. 2009). For the Kuiper belt, we 
assumed an annulus extending from 30 AU to 55 AU, giving a temperature range of 38 - 51 K. 
In both cases, the dust temperatures were calculated assuming blackbody dust grains in radiative 
equilibrium with a Sun-like star. The modeled Ldust/A* values for the zodiacal dust (~ 10 -7 ; 
Nesvorny et al. 2010) and dust in Solar System’s Kuiper belt ( 10 7 ; Vitense et al. 2012) are 

marked with horizontal bars. The latter value is highly dependent on poorly known characteristics 
of the source bodies (Kuiper belt objects); it was revised downward by about a factor of 10 in the 
last two years, based on new knowledge from in-situ detections of dust grains with the Student 
Dust Counter on-board NASA’s New Horizons mission (Poppe et al. 2010; Han et al. 2011). 

For nulling instruments (KIN and LBTI), to convert the null leakage to a relative flux (Td us t/F*), 
one must account for the fact that some dust emission may be removed by the dark fringes, de- 
pending on the exact spatial distribution of the dust. A recent KIN exozodiacal dust survey found 
an average value of 0.4 for the fraction of dust emission transmitted through the fringes, assuming 
unmodified zodiacal- twin disks (Millan-Gabet et al. 2011). The null leakage error must be divided 
by this transmission factor to give the correct uncertainty on the relative flux. However, the cor- 
rect factor assuming a ring-like dust distribution, as we have done here, is not presently known and 
could be as large as 1. Therefore, we have chosen to ignore the transmission factor in calculating 
the sensitivity curves for nulling instruments. 

Over the last several years, high precision visibility measurements with near-IR interferome- 
ters (VLTI/VINCI, CHARA/FLUOR) have been used to detect circumstellar dust around a few 
nearby stars (e.g. di Folco et al. 2007). These observations can detect large amounts of hot dust 
(CHARA/FLUOR 3a sensitivity ~ 5000x a Solar System-twin disk for ~ 1700 K dust; O. Absil, 
personal communication). While these observations are most sensitive to dust interior to habitable 
zones, they raise interesting questions about the origin of large amounts of dust so close to the 
central stars; one possibility is star-grazing comets (Absil et al. 2006). Another new instrument 
with sensitivity to hot dust is the Palomar Fiber Nuller (PFN; Hanot et al. 2011). Recent PFN 
observations at 2.2 /xm of a debris disk system with a hot dust detection from CHARA (Vega) limit 
the dust to < 0.2 AU from the central star (or > 2 AU if there is a significant scattered light con- 
tribution; Mennesson et al. 2011). Unfortunately, we cannot accurately estimate the sensitivity of 
these near-IR instruments with Equation 4, since the assumption that the stellar spectrum is well- 
described by a Rayleigh- Jeans law breaks down; therefore, we have not included CHARA/FLUOR 
or PFN sensitivity curves in Figure 3. 

Examination of Figure 3 shows that previous and current facilities (Spitzer, WISE, KIN, 
Herschel) are relatively insensitive to dust in the habitable zone of Sun-like stars. To date, the 
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most sensitive survey for habitable zone dust around nearby solar-type stars was performed with 
the KIN (Millan-Gabet et al. 2011). While the mean 3 a detection limit for an individual star in this 
survey was 480 x a zodiacal-twin disk (shown in Figure 3), the ensemble of non-detections was used 
to infer a 3cr upper limit on the mean exozodi level of 150 x a zodiacal-twin disk. Unfortunately, 
these limiting values are far too high to reliably estimate the performance of a future mission aimed 
at direct observations of Earth-analog planets. The sole near-term facility sensitive enough to assess 
the habitable zone dust around solar-type stars is LBTI, which should have a 3<r detection limit 
of about 10 x a Solar System-twin disk for ~ 300 K dust. An LBTI survey for exozodiacal dust 
around a set of about 60 - 100 nearby stars is expected to begin in late 2012 3 . 

Turning to dust structures, observations of dust thermal emission are typically unresolved; at 
best, SED modeling provides information only on the radial distribution of the dust. Dust clumps 
caused by an Earth-mass planet orbiting 1 AU from a Sun-like star may be expected to have sizes 
of about 0.2 AU, corresponding to 0.02 arcsec for a star at 10 pc (Dermott et al. 1994). For each 
facility considered above, the spatial resolution at the reference wavelength is given in Table 1. Of 
the current and near-term facilities, only ALMA will have spatial resolution high enough to resolve 
dust clumps as small as those produced by an Earth-mass planet; unfortunately, ALMA is not 
sensitive to warm dust in habitable zones. However, note that more massive planets can create 
larger clumps (Stark &; Kuchner 2008). 

Looking beyond the timeframe covered above, the next generation of ground-based Extremely 
Large Telescopes (ELTs) should begin operation in the early 2020’s. For these > 25-meter apertures, 
detecting and spatially resolving thermal emission from habitable zone dust around nearby stars 
should be a tractable problem. Assuming a Solar System-twin disk, the exozodi flux at 10 ^ mr 
relative to the stellar flux is Td ust /-F* ~ 5 x 10” 5 , or about 85 /uJy for a system at 10 pc; this is well 
within the photometric detection limits of planned ELTs (A. Weinberger, personal communication). 
The spatial resolution of a 25-meter aperture at that wavelength (A / D k, 83 milliarcsec) could allow 
us to resolve thermal emission from regions at radii > 0.8 AU for a system at 10 pc, as long as light 
from the central star is sufficiently suppressed and removed at such a small angular separation. 

So instead of an unresolved measurement of the combined light from the star and dust (all 
the instruments discussed above do this, except the infrared interferometers and ALMA in the 
extended array mode) , an ELT with advanced high-contrast imaging instrumentation could directly 
image the thermal emission from habitable zone dust around nearby stars. In the next few years, 
new coronagraphs operating at near-IR wavelengths on ground-based telescopes should reach 10~ 6 
contrast at small angular separations, using extreme adaptive optics and reference PSF subtraction 
(fuller discussion of coronagraphic imaging appears in the next section). This contrast will also be 
possible at 10 /rm in the next decade; in fact, it will probably be easier to achieve, thanks to the 
higher Strehl ratio possible at the longer wavelength. A simulation of coronagraphic imaging of 
habitable zone dust around a nearby solar-type star with an ELT appears in Figure 4. It seems 


3 http: //nexsci . caltech. edu/missions/LBTI/cfp_keysci . shtml 
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likely that by sometime in the 2020s, ELTs will measure low levels of dust in the habitable zones 
of stars within ~ 10 pc. 


4. High-Contrast Imaging in Scattered Light 

Debris disks may also be detected through coronagraphic imaging of scattered light; such 
observations have high spatial resolution and are best for revealing dust structures (see Figure 5). 
However, current coronagraphy is far less sensitive at all distances from the central stars than 
observations of thermal emission, and cannot image the habitable zones at all. This is due to 
the difficulty of suppressing diffracted and scattered light from the bright central star (using a 
coronagraph) , and then removing residual starlight through subtraction of a reference PSF. The 
accuracy with which this can be done largely depends on the stability of the telescope PSF, whether 
one is using a space-based or ground-based facility. 

PSF subtraction can be improved using differential imaging techniques, which distinguish 
between PSF artifacts and real sources. Some examples are angular differential imaging (ADI), 
which images the target at many rotation angles (e.g. Marois et al. 2006); chromatic differential 
imaging (CDI), which uses the known wavelength dependence of PSF artifacts to identify and 
remove them (e.g. Crepp et al. 2011); and polarimetric differential imaging (PDI), which compares 
images of the target in different polarizations to remove unpolarized starlight (e.g. Quanz et al. 
2011). Note that ADI and CDI are best at searching for faint point-sources and work less well 
on extended structures like debris disks. There are also post-observation techniques that can help 
with PSF subtraction, like locally optimized combination of images (LOCI; Lafreniere et al. 2007). 

To date, the most sensitive scattered light images of debris disks were obtained with HST, 
due to its excellent PSF stability compared to ground-based telescopes, even those with adaptive 
optics (AO) systems. Several HST instruments have had high-contrast imaging capability; they are 
listed in Table 2. This table also lists the Ld us t/L* value for the faintest disk successfully imaged 
in scattered light with each HST instrument. Examination of these values quickly demonstrates 
the much greater sensitivity of thermal emission observations (compare to Figure 3). 

A new generation of coronagraphs behind extreme AO systems on 8-meter-class ground-based 
telescopes is in the development and commissioning phase. These include the High Contrast In- 
strument for the Subaru Next Generation Adaptive Optics (HiCIAO) instrument on the Subaru 
Telescope, the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) instrument on 
the VLT, and the Gemini Planet Imager (GPI) on the Gemini South Telescope. HiCIAO has be- 
gun science operations (e.g. Thalmann et al. 2011), while we expect the other two instruments to 
be commissioned in the next year. Further in the future, all three imaging instruments planned 
for JWST (NIRCam, NIRISS, and MIRI) include high-contrast capabilities; NIRCam and NIRISS 
(formerly known as TFI) will operate in the visible/near-IR and be sensitive to scattered light. 
More information about all these instruments and others may be found in Beichman et al. (2010). 
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Due to the inherent difficulties in assessing instrument performance before commissioning 
and also to the variety of observing and post-observation processing techniques available, it is 
difficult to predict the ultimate performance of these new instruments for high-contrast imaging 
of extended structures like debris disks. A preliminary measurement of the point-source contrast 
achievable using HiCIAO with ADI is available (see Table 2); so far, the faintest disk successfully 
imaged is HR 4796A, one of the very brightest debris disks (shown in the right-hand panel of 
Figure 5). Table 2 also includes predictions for the contrast that will be achievable with GPI and 
JWST/NIRCam; however, these assumed ADI would be used and may not accurately predict the 
instrument contrast for extended sources. A recent prediction of the GPI performance for debris 
disks suggests that compact (1" diameter) disks with L^ust/^* > 2 x 10~ 5 can be imaged, using PDI 
and assuming moderately polarizing grains (B. Macintosh, personal communication). Unlike the 
other instruments, JWST/NIRISS uses sparse-aperture interferometric imaging to perform high- 
contrast imaging; a prediction for its point-source contrast appears in Table 2, but the expected 
extended-source performance is not yet available. 


5. Summary 

1. Exozodiacal dust affects future efforts to directly observe Earth-like planets in the habitable 
zones of nearby stars in two ways: 1) background flux leading to increased noise and 2) dust 
structures causing confusion with unresolved exoplanets. The impact of second problem has 
not been thoroughly studied. 

2. Assuming uniform exozodi surface brightness, a 4-meter telescope aperture, and optical ob- 
serving wavelengths, the counts from the local zodiacal dust plus a Solar System-twin disk 
of exozodi are about 5 times greater than the counts from the Earth observed at quadrature 
orbiting 1 AU from a Sun-like star at 10 pc. 

3. LBTI is the only previous, current, or near-term facility sensitive enough to detect exozodiacal 
dust in the habitable zones of nearby solar-type stars at levels approaching the Solar System 
zodiacal dust level. This facility should begin surveying nearby stars for warm dust in 2012. 
A decade or so in the future, ground-based ELTs may be able to image low levels of habitable 
zone dust around stars within ~ 10 pc. 

4. Dust structures like clumps located far from the central stars may currently be detected with 
high-contrast imaging of light scattered from dust disks. These observations are currently far 
less sensitive than observations of unresolved thermal emission. Cold clumps at large distances 
may soon also be imaged with ALMA. Unfortunately, dust structures in habitable zones are 
likely to prove elusive, although there is a chance that a new generation of ground-based 
coronagraphs may provide some information for the nearest stars. 


This work was performed by members of the Debris Disks & Exozodiacal Dust Study Analysis 
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Group (SAG # 1 ), part of NASA’s Exoplanet Exploration Program Analysis Group (ExoPAG). 
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Table 1. Detection of dust thermal emission: data for various telescope facilities 


Telescope / 
Instrument 

Operation 

Dates 

^obs. 

(nm) 

7 -1 a 

± peak 

(K) 

Resolution b 

(") 

Uncertainty c 

(lcr) 

°>dust/D d 

Ref. 

KIN 

2005 - 

- 12 

8.5 

432 

0.01 e 

0.3% leak error f 

0.003 

1 

Spitzer/IRS 

2003 - 

09 

10 

367 

2.4 

1.5% of star flux 

0.015 

2 

LBTI 

2012 - 


10 

367 

0.05 e 

0.01% leak error f 

0.0001 

3 

WISE/W4 

2009 - 

- 11 

22.1 

166 

12 

3% of star flux g 

0.03 

4 

JWST/MIRI 

2018 - 

- 23 

25.5 

144 

0.9 

2% of star flux 

0.02 

5 

Spitzer/MIPS 

2003 - 

- 09 

70 

52 

18.0 

15% of star flux 

0.15 

6 

Herschel/PACS 

2009 - 

- 13 

70 

52 

5.2 

1.61 mJy 

0.04 h 

7 




100 

37 

7.7 

1.90 mJy 

0.1 h 

7 




160 

23 

12.0 

3.61 mJy 

0.4 h 

7 

ALMA 1 

2012 - 


1250 

3 

0.02 

0.1 mJy 

0.7 h 

8 


a Temperature for which the dust blackbody emission peaks at the observation frequency. 
Instrument is most sensitive to dust near this temperature. 

b FWHM of the instrument PSF at the observation wavelength. 

c Photometric uncertainty at the observation wavelength. 

d Photometric uncertainty relative to the stellar flux at the observation wavelength. 

e l/2 the fringe spacing. Spatial resolution only in the direction perpendicular to the fringe 
pattern, so not applicable for dust clumps. 

f Null leakage error. 

Preliminary value, applicable only for the best cases (D. Padgett, personal communication). 

h Calculated assuming the photometric uncertainty in the previous column and F* at the 
observation wavelength for a Sun-like star at 10 pc. 

Mhe ALMA spatial resolution is for the extended array mode, while the photometric uncer- 
tainty assumes an unresolved observation with the compact array. 


References. [1] Millan-Gabet et al. (2011); [2] Beichman et al. (2006); 

[3] Hinz et al. (2008); [4] User’s Guide to the WISE Preliminary Data Release, 
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http : / /wise2 . ipac . caltech . edu/docs/release/prelim/expsup/wise_prelrel_toc . html; 

[5] Estimated final absolute photometric accuracy (C. Chen, personal communication); 

[6] Bryden et al. (2006); [7] Sensitivity of PACS scan map with on-source integration 

time ~ 360 sec, calculated with HSPOT v6.0.1; [8] Expected full array sensitivity in 

Band 6 with integration time of 60 sec, ALMA Early Science Primer, v2.2 (May 2011), 

http : //almatelescope . ca/ALMA-ESPrimer . pdf. 
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Table 2. High-contrast optical/near-infrared imaging of dust scattered light: instrument 

performance 


Facility / 

Operation 

IWA a Contrast 15 

Faintest Disk Imaged Refs. 

Instrument 

Dates 

(") at 1" 

ID (Ldust/A*) 


HST/STIS 

1997 - 2004, 2009 - 

0.5 

HST/NICMOS 

1997 - 99, 2002 - 08 

0.5 

HST/ACS 

2002 - 07 

1 

Subaru/HiCIAO 

2010 - 

0.15 

Gemini S/GPI 

2012 - 

0.08 

JWST/NIRCam 

2018 - 23 

0.3 

JWST/NIRISS 

2018 - 23 

0.1 


3 x IO -3 

HD202628 

1 x 10 -4 

1, 2 

10 -5 

HD181327 

2 x 10“ 3 

3,4 

10 -5 

Fomalhaut 

8 x 10 -5 

5, 6 

io ~ 4 - 8 c 

HR4796A 

5 x 10 -3 

7, 8 

, io -6 to ~ 7 c 



9 

~ 10“ 5 c 



10 

, io -4 to _5 c 



11 


a Inner working angle (smallest achievable). 

b Relative to peak of unobscured PSF, with reference PSF subtracted. 
c Assuming a point-source. Will probably be worse for extended sources like disks. 

References. — [1] STIS Instrument Handbook, vlO.O, 

http : //www . stsci . edu/hst/stis/documents/handbooks/ currentIHB/cover . html; 

[2] Krist et al. (2012); [3] Schneider & Hines (2007); [4] Schneider et al. (2006); [5] ACS Instrument 
Handbook, vlO.O, http : //www . stsci . edu/hst/acs/documents/handbooks/cyclel9/cover .html; 
[6] Kalas et al. (2005); [7] Suzuki et al. (2010); [8] Thalmann et al. (2011); [9] GPI web page, 
http://planetimager.org/pages/gpi_tech_contrast.html; [10] Krist et al. (2007), [11] STScI 
NIRISS web page, http://www.stsci.edu/jwst/instruments/niriss/ObservationModes/saii 
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Fig. 1. A composite image of (3 Pictoris system showing the edge-on debris disk and the exoplanet 
/3 Pic b (Lagrange et al. 2010). The disk image was obtained with the ESO 3.6-rn Telescope; the 
innermost portion of the disk is obscured. The planet, imaged at two epochs with ESO’s VLT, 
is likely responsible for a warp seen in the inner disk (Heap et al. 2000), revealed as an inclined 
sub-disk in HST ACS coronagraphic images (Golimowski et al. 2006). Image credit: ESO/A.- 
M. Lagrange. 
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Fig. 2. — Impact of exozodiacal dust emission on direct observations of exoEarths. For both plots, 
we assumed the exozodi has a uniform surface brightness. Both x-axes give the exozodi surface 
brightness at the EEID relative to the brightness of a Solar System-twin disk (e in Equations 2 & 
3). Left: Dust counts and imaging exposure time versus exozodi brightness. The left y-axis is the 
background counts from zodiacal plus exozodiacal dust divided by the counts from an Earth-twin 
around a Sun- twin at 10 pc, both calculated for V band (solid lines). The right y-axis is the exposure 
time (in arbitrary units) required to detect the Earth- twin with some S/N, assuming background- 
limited imaging (dashed lines). The counts ratios and exposure times were calculated for three 
different telescope aperture diameters: 2 meters (black lines), 4 meters (red lines), and 8 meters 
(orange lines). Smaller apertures are more sensitive to background emission. Right: The combined 
effect of exozodi emission and r/E art h on the yield of a direct imaging/spectroscopy exoplanet mission. 
The y-axis is the fraction of stars with an Earth-size planet in the habitable zone. The curves were 
created using a simple mission planning code that chooses real stars within 30 pc for observation 
until 5 mission years is reached, assuming a 4-meter telescope, some value of rj Earths and that all 
the stars have the same exozodi level. At each value of 77 Earth > the “Tolerable Exozodi” (x-axis) 
is the largest exozodi level for which the desired mission yield (expected number of exoEarths 
characterized) is achieved. The analysis was performed for three values of mission yield: high 
(dotted line), medium (dashed line), and low (solid line). The smaller 77 Earth is, the lower the 
exozodi level that can be tolerated while still characterizing the desired number of exoEarths. A 
full description of the code appears in Turnbull et al. (2012); however, the general behavior shown 
here should be similar for all direct observation exoEarth missions. 
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Fig. 3. — Sensitivity limits for detection of debris dust around nearby Sun-like stars, for various 
recent (Spitzer, WISE), current (KIN, Herschel) and near-term facilities (LBTI, JWST, ALMA). 
The curves show 3cr detection limits in terms of the fractional dust luminosity (Tdust/-^*) versus 
its temperature. Recent and current facilities are plotted with solid lines, near-term ones with 
dashed lines. The instrument data assumed for these curves appear in Table 1. The temperature 
ranges for dust in two zones around the Sun are shown with vertical bars, calculated assuming 
blackbody grains. The habitable zone (0.8 AU -1.8 AU) is shown in pink, the Kuiper belt (30 AU 
- 55 AU) is shown in light blue. The modeled Tdust /T* values for the Solar System’s Kuiper belt 
dust (~ 10~ 7 ; Vitense et al. 2012) and the zodiacal dust (~ 10 -7 ; Nesvorny et al. 2010) are marked 
with horizontal light blue and pink bars. 
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Fig. 4. — Simulated high-contrast imaging of exozodi thermal emission with a ground-based ELT. 
Left: A model for a face-on debris disk around r Ceti, a Sun-like star at 3.7 pc (Stark & Kuchner 
2008). The disk has a dust abundance equal to that of a zodiacal-twin disk and an embedded 
1 Earth-mass planet orbiting at 1 AU (planet located at 12 o’clock). Image credit: C. Stark (model 
available for download at http://asd.gsfc.nasa.gov/Christopher.Stark/catalog.php). Middle: A 
simulated 10 ^urn imaging observation of the model using an idealized coronagraph designed to 
produce high Strehl ratios (~ 98% at 10 /jrn) on the Giant Magellan Telescope (25-meter aper- 
ture); theoretically, the instrument suppresses the light from the central star by a factor of ~ 10 3 . 
Telescope and sky background emission were included in the simulation. Right: The simulated 
image after subtraction of a reference PSF image, achieving a contrast of ~ 10~ 6 at an angular 
separation of 1 X/D (83 milliarcsec). Slightly different seeing was assumed for the target and PSF 
observations. 
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Fig. 5. — High-contrast imaging of the HR 4796A debris disk in scattered light. For this disk, 
Adust/-^* = 4.7 x 1(T 3 (Moor et al. 2006). Left: Optical wavelength HST/STIS coronagraphic 

image. Image credit: Schneider et al. (2009). Right: Near-IR wavelength Subaru/HiCIAO image. 
Image credit: Thalmann et al. (2011). 




